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We present results for the mass of the rj' meson for two-flavor lattice QCD in the continuum limit, calculated 
on the CP-PACS computer, using an RG-improved gauge action and clover fermion action with tadpole-improved 
Csw. Measurements are made at three couplings corresponding to a « 0.22,0.16,0.11 fm for four quark masses 
corresponding to 171^^/nip « 0.8, 0.75, 0.7, 0.6. The two-loop diagrams are evaluated using a noisy source method. 
Quark smearing for both one- and two- loop diagrams is successfully applied to obtain ground state signals in the 
rj' channel. We obtain m,j/=0.960(87)10'jgg GeV in the continuum limit, where the second error represents the 
systematic uncertainty coming from varying the functional form for chiral and continuum extrapolations. 


1. Introduction 

An accurate lattice calculation of m^i would 
distinguish between QCD and quark models, and 
provide information relevant to the U(l) prob¬ 
lem. Those attempts carried out so far [1], how¬ 
ever, have employed the quenched approximation 
or used only one value of lattice spacing a for full 
QCD. In this work we make a step forward by us¬ 
ing A^f = 2 full QCD at several a, and evaluating 
rriri' in the continuum limit. 

A technical difficulty which is important to 
overcome is the poor signal-to-noise ratio of the 
jf propagator Its error increases quickly 

with the time separation t so that an effective 
mass plateau is not in general observed when 
a point quark source is used. One sometimes 
fits rrijfi from the ratio Grj'{t)/ G,r(i)= 1 — 
Aexp((m,,' —mTr)t), hoping for cancellation of ef¬ 
fects of excited states. However there is no ratio¬ 
nale behind this expectation. We improve upon 
this by employing a smearing method which en¬ 
ables us to obtain ground state signals in the rj' 
channel for all simulation parameters. 

After a description of our method in Sec. 2, 


we evaluate in Sec. 3 m,,' including systematic 
error from chiral and continuum extrapolations. 
A correlation between topology and m,,' which we 
have observed is presented in Sec. 4. 

2. Calculational Method 

We use Af = 2 gauge configurations previously 
generated for our study of the spectrum of or¬ 
dinary hadrons [2], with an RG-improved gluon 
action and clover quark action with tadpole- 
improved Csw We use three couplings [3 = 1.8, 
1.95 and 2.1 corresponding to a ~ 0.22, 0.16, 0.11 
fm, and four hopping parameters k matched to 
correspond to m,r/wp « 0.8, 0.75, 0.7, 0.6. The 
box size LgO « 2.5 fm has also been matched be¬ 
tween /3. Valence and sea quark masses are iden¬ 
tical. We take (uu + dd)/'/2 as an operator for 
the 77 '. 

G^'(t) = Gi(t) — G 2 (t) is calculated on 400- 
800 stored configurations for each (/3, k) pair, 
where Gi{t) = G-K{t) (G 2 (t)) is the contribu¬ 
tion of the one- (two-) loop diagram. We use 
an exponential-like smearing function given in 
Ref. [2] for both Gi and G 2 . One or two quark 
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Figure 1. Relative error of versus -/Vnoise- x- 
scale is proportional to 1/V-^noise- 
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Figure 2. Comparison of rj' effective masses ob¬ 
tained from local and smeared sources at /3 = 2.1 
and mTr/wp « 0.7. niesit) comes from G{t — 
and G{t + i). 

sources are smeared, while sink is always local. 
Errors are estimated by the jackknife method us¬ 
ing bins of 50 trajectories. The lattice scale is 
hxed from Wp. 

After fixing configurations to the Coulomb 
gauge, we estimate G 2 {t) using a noisy source 
method with C/(l) random fields. In order to en¬ 
hance signals, inversions for the quark propagator 
are carried out independently for each (color,spin) 
index. Fig. 1 obtained for the coarsest lattice 
at /3 = 1.8 shows that increasing the number of 
noisy samples from A^noise = 3 to 20 gives no sub¬ 
stantial gain in precision in rup'. We thus choose 
A^noise = 3 for finer lattices at /3 = 1.95 and 2.1. 
At /3 = 1.8, we use results with A^noise = 8, giving 
the smallest error in the chiral limit. 

Fig. 2 compares effective masses of rrip' ob¬ 
tained from one-quark smeared and local sources. 



p=1.8 P=1.95 P=2.1 


Figure 3. Comparison of rup' from direct fit and 
ratio method. 
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Figure 4. Chiral extrapolations at /? = 2.1. 

(Local results plotted are calculated on the fly 
by the volume source method without gauge fix¬ 
ing, and have smaller errors than those calculated 
from noisy sources; these have already been re¬ 
ported in Ref. [3].) Smearing has the clear effect 
of ensuring a plateau from t « 1.5 for all parame¬ 
ter values. Therefore rup/ is determined by fitting 
smeared Gp' for tmin > 1 using a single hyper¬ 
bolic cosine function. We note for this smeared 
data that fitting the ratio Gp' /G^ gives a result 
consistent with the direct fit as shown in Fig. 3. 
The deviation remains within 8%. 

3. rj' Meson Mass in the Continuum Limit 

We test two functional forms for chiral extrap¬ 
olation, 1) {amn'Y = A + R(am,r)^ (Goldstone 
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Figure 5. to,,' vs. a with linear, quadratic and 
constant continuum extrapolations. 


boson type: ‘GB’) and 2) anij^' = A + B^arriT^)"^ 
(non-GB type) and find that both reproduce data 
equally well (see Fig. 4 for example) with compa¬ 
rable x^/df = 0.9—1.1 (GB) and 0.4—2.2 (non- 
GB). GB fit is used to determine central values, 
while non-GB fit is used for error estimation. 

Fig. 5 shows TO,,' as a function of a. Continuum 
extrapolation is carried out linear in a, since 0(a) 
error remains in our action combination. We also 
try a constant plus quadratic form, since most 
0(a) effects may be removed by use of the im¬ 
proved action. Finally, since to,,' hardly changes 
over the finest two lattices, we make a constant 
fit removing the coarsest point. The latter two 
are used to estimate systematic error. 

Systematic errors arising from chiral and con¬ 
tinuum extrapolations are added in quadrature, 
separately for upper and lower sides. The final 
result reads 

TO,,' = 0.960(87)t;;;^^^ GeV (1) 


4. Correlation with Topological Charge 

The large rj' mass compared to the pion octet 
may originate from the topological structure of 
QCD. Using results for topological charge [4] 
Qtopo, we split the measurements of the rj' prop¬ 
agator into two equally sized bins of high and low 
IQtopol, and measure to,,' on each bin. As shown 
in Fig. 6, configurations with high |Qtopo| give 
large values of to,,' for lighter quarks, though in¬ 
terpretation of results remains as open problem. 
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Figure 6. to,,' from configurations with high and 
low IQtopol normalized by the value for all config¬ 
urations. 


5. Conclusions 

Our value to,,' for Ni = 2 full QCD in the con¬ 
tinuum limit, Eq. 1, turns out to be consistent 
with the experimental value, though there still 
remains considerable numerical uncertainty, par¬ 
ticularly on the lower side (—30%), coming from 
the continuum extrapolation. Additionally, our 
result contains systematic error from quenching 
effects of the strange quark and from neglect¬ 
ing mixing with the ss state. Calculations are 
currently being performed analysing mixing be¬ 
tween quark-based states and mass eigenstates for 
N{ = 2 QCD. 
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